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Abstract 

The doping, temperature and energy dependence of the dynamical spin struc- 
ture factors of the underdoped lanthanum cuprates in the normal state is 
studied within the t-J model using the fermion-spin transformation technique. 
Incommensurate peaks are found at [(l±£)7r, 7t], [tt, (1 ± <5)vr] at relatively low 
temperatures with 5 linearly increasing with doping at the beginning and then 
saturating at higher dopings. These peaks broaden and weaken in amplitude 
with temperature and energy, in good agreement with experiments. The the- 
ory also predicts a rotation of these peaks by 7r/4 at even higher temperatures, 

being shifted to [(1 ± 5/V2)ir, (1 ± 8/y/2)ir\. 
74.25.Ha, 74.72.-h, 74.20.Mn 
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In spite of the tremendous efforts dedicated to the studies of anomalous properties of 
high T c superconductors, many important problems still remain open. Among others, the de- 
struction of antiferromagnetic long range order (AFLRO) and appearance of incommensurate 
antiferromagnetism (IAF) in doped cuprates is one of the challenging issues for the theory of 
strongly correlated electron systems. Moreover, the interplay of AF and superconductivity 
in these compounds is of fundamental importance for the high T c theory. Experimentally, 
by virtue of systematic studies using NMR and /iSR techniques, particularly the inelastic 
neutron scattering, rather detailed information on dynamical magnetic properties has be- 
come available now, awaiting an adequate theoretical interpretation. It has been established 
that beyond certain critical doping (~ 3%) the commensurate AFLRO disappears, being re- 
placed by IAF, characterized by incommensurability parameters 5, i.e., the AF Bragg peaks 
are shifted from [11,71} to four points [tt(1 ± 5),ir], [ti, (1 ± 5)n] |L]. For very low dopings S 
varies almost linearly with concentration x, but saturates at higher dopings. These peaks 
broaden and weaken in amplitude as the temperature and energy increase. These features 
are fully confirmed by the data on lanthanum cuprates and have also been found re- 

cently on yttrium cuprates Theoretically there is a general consensus that IAF emerges 
due to doped charge carriers. Several attempts have been made to make this argument more 
precise, including the hole induced frustration ||, stripe formation @, spiral phase f| and 
Fermi surface nesting ||. Based on the phenomenological ansatz of marginal Fermi liquid 
behavior |1(J and tight binding calculation, a detailed fitting of the experimental data was 



attempted ||11|| . Recently experiments show stronger singularities of AF fluctuations (12 



than what is anticipated from the phenomenological models. The proximity to a quantum 
critical point [13| was proposed as an alternative explanation ||12|| . However, to the best of 



our knowledge, no systematic calculations have been performed within the standard strong 
correlation models for the dynamical spin structure factors (DSSF) to confront the exper- 
imental data. Exact diagonalization is limited by system sizes, while the quantum Monte 



Carlo technique faces the negative sign problem for lower temperatures [3. Thus it is 
rather difficult to obtain conclusive results. 



In this paper, using the fermion-spin theory |15| which implements properly the local 



single occupancy constraint, we calculate explicitly DSSF for cuprates within the t-J model 
and reproduce all main features found in experiments including peak position as 

well as temperature and energy dependence. Apart from the ratio t/J (taken to be 2.5), 
there are no other adjustable parameters in the calculations. Moreover, the theory predicts 
the magnetic peaks will be rotated by 7r/4 at even higher temperatures, i.e., being shifted 
to [(1 ± 5/\/2)7r, (1 ± 5/y/2)ir]. To avoid complications due to bilayers we will focus on the 
normal state IAF in lanthanum cuprates. 

We start from the t-J model on a square lattice, 

H = -t £ ClC m<J + h.c. - im £ ClQ a + J^Sr S l+f)} (1) 

if/a i(r if) 

with the local constraint S CT C| .Cj (T < 1, where i) = ±x, ±y, and Sj = CjaCi/2 are spin 
operators with a = (<j x , <y y ,a z ) as Pauli matrices. The single occupancy local constraint 
can be treated properly in analytical form within the fermion-spin theory |T5| based on the 



slave particle approach = h\Sj~, = h\Sf , where the spinless fermion operator hi 

describes the charge (holon) degrees of freedom, while the pseudospin operator Si describes 
the spin (spinon) degrees of freedom. In this representation, the low- energy Hamiltonian of 



the t-J model (1) can be rewritten as JT5| 



H - t^2hl +fj hi(SfS i+fj + Si S+ +fj ) + h\hi + J eff ^(S, • S i+ f,), (2) 



IT] I IT] 



with J eff = J[(l — x) 2 — 2 ], where x is the hole doping concentration, the holon hopping 
parameter = {h\hi +ri ), and and S^ as the pseudospin raising and lowering operators, 



respectively. It has been shown [[HJ that the constrained electron operator can be mapped 
exactly using the fermion-spin transformation defined with an additional projection operator. 
However, this projection operator is cumbersome to handle in the actual calculations, and we 



have not presented it explicitly in Eq. (2). It has also been shown |15|] that such treatment 
leads to errors of the order x in counting the number of spin states, which is negligible 
for small dopings. Within this framework the spin fluctuations only couple to spinons, but 



the strong correlation between holons and spinons is included self-consistently through the 
holon's parameters entering the spinon propagator. Therefore both spinons and holons are 
responsible for the spin dynamics. The universal behavior of the momentum-integrated 
DSSF and susceptibility in the underdoped regime has been calculated within the fermion- 
spin theory |T7[ by considering spinon fluctuations around the mean-field (MF) solution, 



where the spinon part is treated by the loop expansion to the second order. Following Ref. 



17|| , we obtain DSSF as, 

S(k, to) = -2[1 + n B {uj))ImD(k, uj) 



-2[1 + n B (uj)} ; 75 m > ( 3 ) 



gjlmS( 2 )(k,^j 
[uj 2 -uj 2 - B k Re^\k,uj)} 2 + [B k lm^>(k,uj)} 2 ' 

where the full spinon Green's function, Z} -1 (k, uj) = D^~ 1 (k, uj) — E^(k, uj), with the 
MF spinon Green's function @, D^~ 1 (k,u) = (uj 2 -uo 2 k )/B k , while ImE( 2 ) (k, uj) and 



Re£^ 2 )(k, uj) are the imaginary and real parts of the second order spinon self-energy, respec- 
tively, obtained from the holon bubble, 



E^(k, u) = - — ^(7 p / +p+fc + 7fe _/ ( 



N J , 2uJ k+p \UJ + t; p+p i — £ p i + UJ k+p 

F 2 (k,p,p') \ 



(4) 



u + £p+p' — ip 1 ~ ^k+p J 

where 7k = (1/Z) J2fi e lk '^, Z is the coordination number, B k = A[2xz(£Jk — 1) + x(lk — e)], 
A = 2ZJ cS , e = 1 + 2t0/J eff , Fx(k,p,p') = n F (£ p+p >)[l - n F (£ p >)) + [1 + n B (uj k+p )][n F (^ p/ ) - 
n F (£ p+p >)], F 2 (k,p,p') = n F (£ p+p >)[l - n F (£ p r)) - n B (uj k+p )[n F (£ p <) - n F {£p+p>)], n F (£ k ) and 
n B (uj k ) are the fermion and boson distribution functions, respectively, the MF holon ex- 
citation £ fc = 2Ztx~ik + A*, and the MF spinon excitation, uj 2 = A 2 (A 1 ^l + A 2 r y k + A 3 ) 
with A x = ae( X /2 + e Xz ), A 2 = e[(l - Z)a(e X /2 + Xz )/Z - a(C z + C/2) - (1 - a)/(2Z% 
A 3 = a(C z + e 2 C/2) + (l — a)(l + e 2 )/(4Z) — ae(x/2 + eXz)/Z, the spinon correlation functions 
X = (S?Sr + ~), Xz = (S!Sf +fj ), C = (1/Z 2 ) E^(5&,S^>, and C z = (1/Z 2 ) E^(Sf +fj S^,). 
In order to satisfy the sum rule for the correlation function (S/~Sf) = 1/2 in the absence of 
AFLRO, a decoupling parameter a has been introduced in the MF calculation, which can 
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be regarded as the vertex correction |l^,|T£|. These MF order parameters Xi Xz, C z , 0, 
and decoupling parameter a have been determined by the self-consistent equations. 

Of course, at vanishing dopings the AFLRO gives rise to a commensurate peak at 
[1/2, 1/2] (hereafter we use the units of [2n, 2tt]), which is not presented here for the sake of 
space. Instead, we plot DSSF S(k,u) in the (k x , k y ) plane at doping x = 0.06, temperature 
T = 0.1J and energy u = 0.05 J for t/J = 2.5 in Fig. 1. The commensurate peak is split into 
four IAF peaks at [(l±5)/2, 1/2] and [1/2, (l±5)/2\. The calculated DSSF S(k,w) has been 
used to extract the doping dependence of the incommensurability parameter S(x), defined 
as the deviation of the peak position from the AF wave vector [1/2, 1/2], and the result 
is shown in Fig. 2 in comparison with the experimental data |4| taken on L^-zSr^CuC^ 
(inset). 5(x) increases almost linearly with the hole concentration in the low-doping regime, 
but it saturates at higher dopings, in full agreement with experimental data. 

For a better understanding of the IAF we have made a series of scans for S(k, uj) at 
different temperatures and energies, and the result for doping x = 0.06, t/J = 2.5 at 
T = 0.1J and uj = 0.1J is shown in Fig. 3. Comparing it with Fig. 1 for the same 
set of parameters except for uj = 0.05J, we see that at relatively low temperatures (T = 
0.1J), although the positions of IAF peaks are almost energy independent, these peaks are 
broadened and suppressed with increasing energy, and tend to vanish at high energies. This 
reflects that the spin excitations are rather sharp in momentum space at low temperatures 
and energies, then the linewidth, or the inverse lifetime increases with increasing energy, in 
full agreement with experiments pj,|12||. Now we turn to discuss the temperature dependence 
of S(k,uj). S(k,u>) at x = 0.06 for t/J = 2.5 and uj = 0.05 J at temperature T = 0.5 J is 
plotted in Fig. 4. To our big surprise, comparing it with Fig. 1 for the same set of 
parameters except for T = 0.1J, we find that, apart from the suppression of the peak weight 
with temperature as anticipated, the positions of IAF peaks are temperature dependent, i.e., 
these peaks deviate from [(1 ± S)/2, 1/2] and [1/2, (1 ± 5)/2] with increasing temperature, 
and are rotated by n/4 in the reciprocal space about (1/2,1/2) at higher temperatures 
(T > 0.5 J), being shifted to [(1 ± 6/y/2)/2, (1 ± 8/y/2)j2\. Up to now most experimental 



data show that the positions of IAF peaks in L^-xSr^CuC^ and L^SrCuC^+z |20 

are located at [(1 ± 6)/2, 1/2] and [1/2, (1 ± 5)/2] in the underdoped regime, but these data 
in the normal-state are obtained at relatively low temperatures (near the superconducting 
transition). However, a strong temperature dependence S(k,u) in La2- x Sr x Cu04 has been 
observed §0, namely the weight of IAF peaks at [(1 ± 5)/2,l/2] and [1/2,(1 ± 6)/ 2] 
is suppressed severely with increasing temperature, whereas the weight is increasing with 
temperature |L2f at [(l±<5/2)/2, (l±<5/2)/2]. This tendency is consistent with our theoretical 
predictions. Experiments at even higher temperatures are required to check our predictions 
explicitly. 

Now we give some physical interpretation to the above obtained results. As seen from 
Eq. (3), the DSSF has a well-defined resonance character. S(k, u) exhibits a peak when 
the incoming neutron energy u is equal to the renormalized spin excitation E\ = io\ + 
B k Re^\k,E k ), i.e. W(k s ,u) = [uo 2 -uj 2 kg -5 fca Re£( 2 )(k 5 , u)f = (w 2 -£fj 2 ~ for certain 
critical wave vectors (positions of IAF peaks). The height of these peaks is determined by 
the imaginary part of the spinon self-energy, i.e. l/ImE^^k^, to). Near half-filling, the spin 
excitations are centered around the AF wave vector [1/2, 1/2], so the commensurate AF peak 
appears there. Upon doping, the holes disturb the AF background. Within the fermion- 
spin framework, as a result of self-consistent motion of holons and spinons, IAF is developed 
beyond certain critical doping, which means, the low-energy spin excitations drift away from 
the AF wave vector, or the zero of W(k.s,u)) is shifted from [n,ir] to k^. As seen from Eq. 
(3), the physics is dominated by the spinon self-energy renormalization due to holons. In 
this sense, the mobile holes are the key factor leading to IAF. As seen from Fig. 5, function 
W(\t, uj) has a rather deep valley along a circle of radius 5 around [1/2, 1/2]. However, if we 
enlarge the scale very significantly (by a factor of 300!!), as shown in Fig. 6, there is a strong 
angular dependence with actual minima (not exactly zero due to precision limitations) at 
[(1 -5)/2, 1/2] and [1/2, (1 -5)/2] for T = 0.1 J. These are exactly the positions of the IAF 
peaks determined by the dispersion of very well defined renormalized spin excitations. Since 
the height of the IAF peaks is determined by damping, it is fully understandable that they 
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are suppressed as the neutron energy u and temperature are increased. The novel result of 
this paper, namely the shift of the position for IAF peaks as the temperature increases is also 
due to the same reason. To demonstrate this point in Fig. 6, we plot the function W(k, u) 
along the arc from k = [(1 - S)/2, 1/2] via [(1 - 6/y/2)/2, (1 - 6/y/2)/2] to [1/2, (1 - 5)/2] 
at doping x = 0.06 for t/J — 2.5 and uj = 0.05J for different temperatures T = 0.1J (solid 
line), T = 0.2 J (dashed line), T = 0.4 J (dash-dotted line), and T = 0.5 J (dotted line) which 
shows clearly the temperature dependence of the minima of W(\s.,u) - kj. This means the 
self-energy correction due to holon motion is temperature dependent. From the physical 
point of view, this is very reasonable. It would be of great interest to check explicitly this 
prediction in neutron experiments at much higher temperatures. 

To conclude we have shown very clearly in the paper that if the local single occupancy 
constraint is treated properly (as done in the fermion-spin theory) and the strong spinon- 
holon interaction is taken into account, the t-J model per se can correctly reproduce all 
main features of IAF in underdoped cuprates, including the doping dependence of the IAF 
peak position and the energy as well as temperature dependence of the amplitude of these 
peaks, without using adjustable parameters. We believe these are universal features of the 
underdoped cuprates, as shown by experiments on L^-zSr^CuC^ La 2 Cu0 4+x PP] 

and YBa2Cu3C>7_ x . f5j. There might be additional features due to bilayer splitting in the 
band structure [0,0], and related theoretical results will be presented elsewhere |22| . The 
theory also predicts a rotation of IAF peak position at very high temperatures which should 
be verified by future experiments. 

Finally, we would like to mention that the influence of the additional second-neighbor 
hopping t' on the IAF and momentum-integrated dynamical spin susceptibility of the t-J 
model has been discussed within the fermion-spin theory. It has been shown [^3| that for 
small values of t' the qualitative behavior of the IAF and integrated dynamical susceptibility 
of the t-t'-J model is the same as obtained from the present t-J model. 
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FIGURES 

FIG. 1. The dynamical spin structure factor in the (k x ,k y ) plane at doping x = 0.06 at tem- 
perature T = 0.1J and energy oj = 0.05 J for parameter t/J = 2.5. 

FIG. 2. The doping dependence of the incommensurability 5{x) of the antiferromagnetic fluc- 
tuations. Inset: the experimental result on La2_ a; Sr :E Cu04 taken from Ref. [4]. 

FIG. 3. The dynamical spin structure factor in the (k x ,k y ) plane at doping x = 0.06 for 
parameter t/J = 2.5 and energy oj = 0.1J at temperature T = 0.1J. 

FIG. 4. The dynamical spin structure factor in the (k x ,k y ) plane at doping x = 0.06 for 
parameter t/J = 2.5 and energy to = 0.05 J at temperature T = 0.5 J. 

FIG. 5. Function W(k,uj) in the (k x ,k y ) plane at doping x = 0.06 for parameter t/J = 2.5 
and energy w = 0.1 J at temperature T = 0.1 J. 

FIG. 6. Function W(k,w) from k x = [(1 - S)/2, 1/2] via k 2 = [(1 - S/y/2)/2, (1 - S/y/2)/2\ 
to k3 = [1/2,(1 — 5)/2] at doping x = 0.06 for parameter t/J = 2.5 and energy u = 0.05J at 
temperatures T = 0.1J (solid line), T = 0.2J (dashed line), T = 0.4J (dash-dotted line), and 
T = 0.5 J (dotted line). 
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